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ABSTRACT 

We extend the concept of galaxy environment from the local galaxy number density to the grav- 
itational potential and its functions like the shear tensor. For this purpose we examine whether or 
not one can make an accurate estimation of the gravitational potential from an observational sample 
which is finite in volume, biased due to galaxy biasing, and subject to redshift space distortion. Dark 
halos in a ACDM simulation are used in this test. We find that one needs to stay away from the 
sample boundaries by more than 30/i^^Mpc to reduce the error within 20% of the root mean square 
values of the potential or the shear tensor. The error due to the galaxy biasing can be significantly 
reduced by using the galaxy mass density field instead of the galaxy number density field. The error 
caused by the redshift space distortion can be effectively removed by correcting galaxy positions for 
the peculiar velocity effects. We inspect the dependence of dark matter halo properties on four en- 
vironmental parameters; local density, gravitational potential, and the ellipticity and prolateness of 
the shear tensor. We find the local density has the strongest correlation with halo properties. This 
is evidence that the internal physical properties of dark halos are mainly controlled by small-scale 
physics. In high density regions dark halos are on average more massive and spherical, and have 
higher spin parameter and velocity dispersion. In high density regions dark halos are on average more 
massive and spherical, and have higher spin parameter and velocity dispersion. We also study the 
relation between the environmental parameters and the subtypes of dark halos. The spin parameter 
of satellite halos depends only weakly on the local density for all mass ranges studied while that of 
isloated or central halos depends more sensitively on the local density. The gravitational potential and 
the shear tensor have weaker correlations with halo properties, but have environmental information 
independent of the local density. 

Subject headings: galaxies: clusters: general-galaxies: evolution-galaxies: formation- galaxies: fun- 
damental parameters - galaxies: environment- galaxies: statistics 



1. INTRODUCTION 

One of recent developments in the study of galaxy for- 
mation is quantitative understanding of the dependence 
of galaxy properties on environment (Park et al. 2007; 
Hwang & Park 2009 among others) . It has already been 
noticed since the 1930's that galaxy luminosity and mor- 
phology depend on the local density: the high density re- 
gions preferentially harbor more luminous and early mor- 
phological type of galaxies (Hubble & Humason 1931). 
In the beginning of the studies on the environmental ef- 
fects on galaxy properties (on the simulation side, see 
Lemson & Kauffmann 1999; Gao et al. 2005; Jing et 
al. 2007; Hahn et al. 2007; Maccio et al. 2007; on the 
observation side, see Blanton & Berlind 2007; Park et 
al. 2007; Cervantes-Sodi et al. 2008; Skibba et al. 2009; 
Blanton & Moustakas 2009), the environment was distin- 
guished according to the large-scale structure where the 
galaxies under study are located. For example, compar- 
ative studies for galaxies located within massive galaxy 
clusters, groups, or voids were carried out (Oemler 1974; 
Dressier 1980; Postman & Geller 1984; Rojas et al. 2004 
among many others). 

Another trend of the same kind of study used contin- 
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uous parameters that measure the local galaxy number 
density. Various kinds of smoothing kernel were used to 
estimate the local number density from galaxy positions 
in redshift space. The most popular one is the truncated 
cylindrical cone which is motivated from the fact that 
massive clusters appear as Fingers of God in redshift 
space (Hogg et al. 2004; Kauffmann et al. 2004; Kuehn 
& Ryden 2005; Reid & Spergel 2009). The Gaussian and 
spline filters are also often used. These filters can have a 
fixed size or can vary in size to include a fixed number of 
galaxies (Park et al. 1994; Monaco & Efstathiou 1999; 
Park et al. 2007; Kim et al. 2008). 

Even though these two approaches look quite differ- 
ent, they essentially use the local galaxy number density 
to distinguish among different environments. Extension 
of the concept of environment beyond the 'local number 
density' has been started by several authors. Park, Gott 
& Choi (2008) and Park & Choi (2009) used the mass 
density due to the galaxy plus dark halo systems as a new 
environmental parameter. They also divided the 'local' 
density into the large-scale background mass density and 
the small-scale density attributed to the nearest neigh- 
bor galaxy. The mass density is estimated from galaxy 
luminosity and mass-to-light ratios. It turned out that 
the environment set up by the nearest neighbor was criti- 
cally important in determining galaxy properties. Lee et 
al. (2009) used galaxy luminosity density and local color 
(difference between luminosity densities in two bands) as 
environmental parameters. Lee & Lee (2008) inspected 
the relation between the ellipticity of the tidal shear and 
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galaxy morphology. 

It is expected that some galaxy properties depend on 
the local galaxy number/mass density sensitively. Color 
and recent star formation activity may be such prop- 
erties. However, the root-mean-square (RMS) displace- 
ment of mass is about 10/i~^Mpc till the present epoch 
in the ACDM model best fit to the Wilkinson Microwave 
Anisotropy Probe (WMAP) 3-year data (Park & Kim 
2009), and therefore the local density at the present lo- 
cation of galaxies cannot fully represent the environment 
where galaxies formed and evolved. This is true particu- 
larly for intermediate and high density regions since their 
sizes are typically only a few Mpcs. Therefore, it may be 
useful to consider environmental parameters other than 
local density to understand the environmental effects on 
galaxy formation. 

A theoretically motivated environmental parameter is 
the gravitational potential. Dark matter and baryons are 
expected to fall into the deep gravitational potential well 
to form massive objects. Since the gravitational poten- 
tial field picks up the fluctuation power at scales much 
larger than those of the density field, the correlation be- 
tween them will not be perfect. And it will be interesting 
to see how galaxy properties are related with the 'local' 
gravitational potential. 

Furthermore, it is expected that galaxy angular mo- 
mentum is generated from the large-scale gravitational 
shear force. According to tidal torque theory (Hoyle 
1951; Peebles 1969; Doroshkevich 1970; White 1984; Lee 
& Pen 2000; Vitvitska et al. 2002; Porciani et al 2002), 
the origin of galactic angular momentum is originated by 
the tidal torque operating on primordial gas lump that 
will form a galaxy. The torque is given by r = T x /, 
where T is a shear tensor generated by external material 
and / is a moment of inertial tensor of material being 
torqued. Navarro et al. (2004) presented supporting ev- 
idence for the theory with direction of galaxy rotation 
axis. An accurate estimation of the tidal shear tensor 
will enable one to verify if the tidal torque theory is re- 
ally responsible for galaxy spin (Porciani et al. 2002; Lee 
& Pen 2002). 

In this paper we will study how accurately one can es- 
timate the gravitational potential and its functions from 
a simulated sample of galaxies. The error sources in this 
estimation are divided into three categories: 

1. finite volume of the survey, 

2. galaxy biasing, and 

3. redshift space distortion. 

We then inspect the dependence of dark matter halo 
properties on various environmental parameters includ- 
ing the 'local' gravitational potential. It is hoped that 
a generalization of environmental parameter beyond the 
local density allow us to better understand galaxy for- 
mation and evolution. 

2. METHOD 

2.1. Simulation 

We will examine how acciirately one can estimate the 
gravitational potential and its derivatives when an ob- 
servational sample is given. For this purpose we use a 
set of mock galaxies identified as dark matter halos in 
an N-body simulation of the universe. The simulation 
we use here adopted the cosmological parameters mea- 



sured from the WMAP 3-year data (Spergel et al. 2007), 
which are Ha = 0.762, = 0.238, 0;, ^ 0.042, = 
0.958, /i = 0.732, and ag = 0.761, where QA,^m,^b are 
density parameters due to cosmological constant, matter, 
and baryon, respectively, Ug is the slope of the primor- 
dial power spectrum, and as is the RMS fluctuation of 
mass in an 8ft-~^Mpc radius spherical top hat. It ran 
2048^ Cold Dark Matter (CDM) particles whose initial 
conditions are generated on a 2048^ mesh in accordance 
with the ACDM power spectrum. The simulation was 
started at Zi = 47 taking 1880 global time steps till the 
present epoch. The physical size of the simulation cube 
is 1024/i~^Mpc. We used a parallel particle-mesh (PM) 
-|- tree N-body code (Dubinski et al. 2004) to increase 
the spatial dynamic range. The gravitational force soft- 
ening length is set to 50/i~^kpc which is 0.1 times the 
mean separation between CDM particles. The particle 
mass is 9.6 x h-^lO^MQ. 

We identified the gravitationally bound, tidally stable 
dark matter halos (physically self-bound or PSB halos) 
from the CDM particle data at 2 = (Kim & Park 
2006). At the first step the Fricnd-of Friend (FoF) al- 
gorithm is used to search for dark halos adopting the 
connection length of 0.2 times the mean particle separa- 
tion. Then subhalos are identified within each FoF par- 
ticle group taking into account the gravitational binding 
energy with respect to the local density maximum and 
the tidal force from other more massive subhalos if they 
exist. The minimum mass halos contain 30 member par- 
ticles, and the halo mass function is accurate down to 
Mh = 2.9 X 10^^h~^MQ. The mean separation between 
dark halos is 4.3ft-~^Mpc. The resulting dark matter halo 
sample consists of three types: isolated halo which does 
not overlap with other halo (but it can have close neigh- 
bor halos), central halo which is the most massive halo 
in each group of halos, and satellite halo which is not the 
central one in each group of halos. 

Galaxies are identified with the halos. This dark 
matter halo-galaxy one-to-one correspondence model de- 
scribes the observed galaxy distribution quite accurately 
(Kim et al. 2008; Gott et al. 2008; Gott et al. 2009). 
Note that our dark halos are not the commonly used 
FoF halos to which the Halo Occupation Distribution 
prescription is usually applied to statistically distribute 
galaxies (Zheng et al. 2008 and references therein). Our 
dark halos can be central galaxies or satellites, which are 
the direct results of the N-body simulation. To match 
the dark halo sample with observed galaxy samples Kim 
ct al. (2008) adjusted the halo mass threshold, mak- 
ing the mean number densities of halos and galaxies the 
same. The resulting halo sample was to be compared 
with a volume-limited sample of galaxies brighter than 
an absolute magnitude threshold. 

2.2. Potential Calculation 

In the simulation we know the positions of all CDM 

particles within the simulation cube with periodic bound- 
aries. Therefore, the gravitational potential can be cal- 
culated through the Poisson equation 

$k = — p-pomk = -2~fc2~'^™k, (1) 

where p is the mean density, 6 is the overdensity, and we 
decomposed the gravitational potential and the density 
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fields into Fourier modes. This true gravitational poten- 
tial is going to be compared with those obtained by using 
galaxies (i.e. dark halos). 

In practice, we cannot observe dark matter particles, 
but can observe only galaxies. We adopt a simple method 
to obtain the gravitational potential and its derivatives 
from a galaxy redshift sample. Wc assume galaxies are 
locally biased tracers of the underlying mass: Sg{x) = 
bSmix), where b is the bias factor. Then the Poisson 
equation in the Fourier space is 



3 f^-jy;^ ft 
'2 bk^ 



(2) 



The proportionality factor between $k and fc~^^gk is not 
important in our study because each field will be normal- 
ized by its RMS value. In section 3.2 we will compare the 
potential calculated from the galaxy distribution with 
that from the dark matter distribution. 

To use equations (1) and (2) we first calculate the 
density field from matter particle or galaxies. Wc as- 
sign them on a cubical mesh using the Triangular- 
Shaped Cloud algorithm (Hockney & Eastwood 1981), 
and Fourier transform the density array to get Smk or 
(5gk- The gravitational potential in real space is then 
obtained from equation (1) or (2), and is interpolated 
at each location of galaxies. The gradient of $ and the 
shear tensor are calculated by finite differencing $ in real 
space. 

Near a galaxy at x = the gravitational potential can 
be given by the Taylor expansion 



$(x) 



(3) 



where are the components of local acceleration, and Aj 
arc the cigcn values of the shear tensor Tij = didj^. The 
coordinate axes are assumed aligned with the principal 
axes of the shear tensor. The trace of the shear tensor 
at the galaxy is 



V2$ = ^Ai=47rG(p(0)-p). 
Bardeen et al. (1986) introduced the parameters 



Al — A3 

^Ea7 



Al — 2A9 



,p ■ 



2EA, 



(4) 



(5) 



characterizing the ellipticity and prolateness of a (den- 
sity) distribution. Lee & Lee (2008) adopted the e 
parameter in their study on tidal shear dependence of 
galaxy morphology. However, galaxies typically are not 
located at extrema of potential field unlike the den- 
sity field at the smoothing scale we will adopt {Rq ~ 
6/i~^Mpc), and the sum of eigenvalues can be zero, mak- 
ing the above parameters undefined. This is particularly 
true for large-scale density and potential fields where the 
galaxy-scale density peaks are not resolved. To avoid 
this problem we adopt the asymmetry and prolateness 
parameters 



-B = Al — A3, P = Al — 2A2 -|- A3 



(6) 



normalized by their RMS values. These parameters are 

correlated with the local density through Equation (4). 
Dependence of galaxy properties on the shear tensor 
should be then studied at fixed local density. 



2.3. Peculiar velocity correction 

In section 3.3 we will study how the potential and its 

functions are biased because they are estimated from the 
galaxy distribution in redshift space. It is shown that 
most of the redshift space distortion effects can be re- 
moved by making a linear estimation of the peculiar ve- 
locity from the redshift space galaxy distribution. The 
linear regime peculiar velocity is calculated as follows. 
The linearized continuity equation is 



95 
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where v = dx/dt is the peculiar velocity in comoving 
space, 6{x,t) is the density contrast 6 = {p — p)/p, and 
D{t) is the linear growth factor. If we define the velocity 
potential (j) by dx/dD = V(/), we obtain 

V • v = DV^(j). (8) 
Combining the equations (7) and (8) yields 

V^f/. = -S/D. (9) 
The peculiar velocity is then given by 

where a{t) is the expansion parameter, / = 
dlogD/dloga, and H{t) = a/ a is the Hubble param- 
eter. Substituting Equation (9) into Equation (10) in 
Fourier space, one finally obtains the relation between 
density and peculiar velocity 



dx 



(11) 



We actually use a biased tracer of matter field, namely 
galaxies. If the galaxy number density fiuctuation is re- 
lated with the matter fiuctuation by 6g = b6, the factor 
in Equation (11) is replaced by Sg^k/b. 

3. RESULTS 
3.1. Boundary effects 

The gravitational potential would be correctly calcu- 
lated if the mass distribution were given over the infinite 
space. In practice, any observational sample covers a fi- 
nite volume of the universe with complicated boundaries. 
In this section we will study the eff'ects of the missing 
data beyond the sample boundaries on gravitational po- 
tential and its derivatives. 

We first calculate the reference potential from all dark 
halos in the whole simulation cube. Since we are inter- 
ested in large-scale environment, we smooth the density 
field by a 6h~^Mpc Gaussian. This scale is also moti- 
vated by the fact that the mean separation between 
galaxies is about 6ft~^Mpc (cf. the 'Best' sample of Park 
et al. 2005; see also Choi et al. 2007). We then calculate 
a potential using only the dark halos in the half-size cube 
near a corner of the simulation box. The remaining part 
of the simulation cube is replaced by the mean dark halo 
mimber density estimated from the half-size cube. To 
inspect the boundary effects further, we also calculate a 
potential from the dark halos in the quarter-size cube at 
the same corner. 

In Figure 1 we compare the three gravitational po- 
tentials in a slice passing through a ^-coordinate of 
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Fig. 1. — EUipticity of the gravitational shear tensor calculated 
from the whole simulation cube data of size 1024 /i~^Mpc (top), 
from a half-size subcube (middle), and from a quarter-size subcube 
(bottom). Shown are constant shear magnitude contours at four 
threshold levels in half-size slices. 



z — 128ft.^^Mpc of the simulation cube. The shoe in- 
cludes the center of the quarter cube. The top panel 
shows the shear ellipticity field E/a^ in the slice calcu- 
lated from all dark halos in the simulation. We show 
only the bottom left corner of the slice where a compar- 
ison with the half cube result is possible. The contour 
levels correspond to la high (red contours), mean, 1.5a 
low, and 2.5a low (black) shear ellipticity. In the middle 
panel the shear ellipticity field calculated from the dark 
halos in the half cube is shown. The contour levels are 
the same. It is clear that the shear field of the half cube is 
extraordinarily similar with that of the full cube. This is 
particularly true near the center of the half cube, namely 
at the position (256, 256). However, large differences are 
also observed near the boundaries. The bottom panel 
shows the shear field of the quarter cube. It again shows 




50 100 150 200 250 



d (h-iMpc) 

Fig. 2. — (top) The Root-Mean-Square (RMS) error in poten- 
tial estimated from half-size (solid line) and quarter-size (dashed 
line) subcubes as a function of distance from cube surfaces. The 
error is normalized with respect to the RMS variation of potential, 
(middle) Same but for the a;-component of the potential gradient, 
(bottom) Same but for the ellipticity of the shear tensor. 



that the shear ellipticity of the quarter cube agree quite 
well with those of the other ones, particularly near the 
center of the quarter cube. It can be seen that there is a 
weak shear field extending beyond the sample boundary. 
The shear field of the half cube also has such leakage 
outside the boundaries (truncated in the middle panel). 

To quantitatively estimate the accuracy of the poten- 
tial and its derivatives obtained from subcubes we cal- 
culate the RMS differences of the fields within a cubical 
shell centered on the center of the subcubes. The solid 
lines in the top panel of Figure 2 are the difference be- 
tween the potential from the full cube and that of the 
half cube within 2 /i^^Mpc-thick cubical shell centered at 
(x,y,z) = (256/i-iMpc,256/i-iMpc,256/i-iMpc). Size 
of a cubical shell is 2(257/i~^Mpc — d), where d is the dis- 
tance of the shell from the surface of the subcube. For 
example, d — l/i~^Mpc corresponds the outermost shell 
of 512/i~^Mpc size. The dashed line is for the quarter 
cube. The top panel shows how the error in the potential 
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Fig. 3. — Correlation between the potential estimated from the 
dark halo number density ($DHn) '^^<i that from dark matter distri- 
bution (<&dm)- Dark halos with mass larger than 2.9 X 10"/i-1Mq 
are used in the potential calculation. 

varies as one moves into the subcubes. Lines and error 
bars are the mean and RMS differences obtained from 
eight subcube rcstilts. The difference monotonically de- 
creases as the distance from boundaries increases. Since 
potentials can be added by a constant without altering 
physics, the difference of about 0.14cr<i> near the center 
can be ignored in the case of half cube. The figure tells 
that, if the potential field is to be obtained with error 
less than say 20% of its RMS value, one should stay more 
than 54/i~^Mpc from all sample boundaries. The quarter 
cube results show a much larger variance that the half 
cube results. 

A similar trend is seen for the potential gradient as 
shown in the middle panel of Figure 2. The x-component 
of the potential gradient mimics the linear radial peculiar 
velocity field very far from an observer. Due to lack of 
data outside subcubes the error in gradient is large near 
boundaries but monotonically decreases as d increases. 
The difference is smaller if one ignores the large-scale 
gradient on the scale of the subcube. 

The bottom panel of Figure 2 compares among the 
shear ellipticity from the full, half, and quarter cubes. 
It is important to note that the error drops as the dis- 
tance from boundaries increases essentially in the same 
way for the half and quarter cube cases. The error 
reaches 20% of the RMS shear ellipticity at the distance 
of d = 33h~^Mpc. The depth of this buffer region must 
be a function of the amplitude and shape of the power 
spectrum. We conclude that one needs to have an obser- 
vational sample much larger than 60/i~^Mpc to secure 
the regions where the shear field has error less than 20% 
of its RMS fluctuation. 

3.2. Effects of Galaxy Biasing 

In the prcndous section we used the mimber density of 
dark halos to calculate the potential field as it is com- 
mon to use the galaxy number density to define environ- 
ment of galaxies. However, it is the matter field that 
determines the true gravitational potential field, and one 
needs to understand the relation between the potential 
from mass and those from mass tracers. In this section 
we assume that each and every PSB dark halo contains 
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Fig. 4. Correlation between the shear tensor magnitude cal- 
culated from the dark matter distribution (a;-axis) and those from 
dark halo number (top) and dark halo mass (middle and bottom) 
distributions. Dark halos with mass larger than 2.9 X W^^h'^MQ 
are used in the top and middle panels. The mass cut is raised to 
1 X 10^^ h~^MQ in the bottom panel. 
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Fig. 5. — Correlation coefficient between the shear tensor magni- 
tude calculated from the dark matter distribution and those from 
dark halo number (dashed line) and dark halo mass (solid line) 
distributions as a function of halo mass cut. 
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one galaxy, and compare the potential calculated from 
their number density or halo mass density with that from 
CDM particles in the simulation. 

In Figure 3 we compare the potential from galaxy num- 
ber density (y-axis) with the correct potential from CDM 
particles. Dark halos more massive than 2.9 x 10^^ Mq/H 
are used for the number density calculation. Their mean 
separation is 4.6ft.^^Mpc, which is equal to that of the 
SDSS galaxies brighter than Mr = -19.5 + 5\ogh (Choi 
et al. 2007). The correlation coefficient between the two 
potential fields is fairly high (r = 0.988). At $dm = 
the potential from the galaxy number density has a dis- 
persion of 0.20 times the RMS potential value. 

Figure 4 compares between the magnitudes of the 
traceless shear tensor Tij = didj^ — -1(5^^X7^$, estimated 
from matter density, halo number density, and halo mass 
density fields. Top panel shows that the correspondence 
between the shear fields from the galaxy number density 
and matter density is not so good. The correlation coeffi- 
cient is only 0.858, and the dispersion of the shear magni- 
tude obtained from galaxy number density is 0.47 times 

the RMS shear magnitude at \T\dmI ^TS^dmY''^ = 2- 

Accuracy in potential and shear can be greatly im- 
proved if the halo mass is used to weigh galaxies and the 

halo mass density field, instead of the number density 
field, is used to calculate the potential. The second and 
third panels of Figure 4 demonstrate such improvement 
when the halo mass threshold is set to 2.9 x 10^^ Mq/Zi 
and 1.0 X IO^^Mq/Zi, respectively. The latter objects 
corresponds to the SDSS galaxies brighter than Mr ~ 
—20.4 + 51og/i, close to that of the M^ galaxies (Choi 
et al. 2007), in the sense that their mean separations 
(6.6/i^^Mpc) are the same. In the middle panel, at 

\T\dm/{\T\dm?''^ = 2, the dispersion in \T\dh is only 
0.17 times the RMS shear magnitude, an improvement 
by almost a factor 3. The correlation drops as the halo 
mass threshold increases, but it still remains quite good 
for massive dark halos with M > 1.0 x lO^^Mo/Zi. 

Figure 5 shows how the correlation changes as a func- 
tion of the halo mass threshold. The solid line is the 
case when the halo mass density field is used, and the 
dashed line is for the halo number density. The cor- 
relation drops rapidly as the threshold increases above 
10^^/i~^M©. We conclude that the gravitational poten- 
tial can be estimated quite accurately by using the ob- 
served galaxy distrilnition, but accuracy can be greatly 
improved if the total mass associated with the dark halo 
plus galaxy system is used to weigh galaxies in the shear 
tensor calculation. 

When one chooses to use the mass field instead of the 
number density field, one needs to adopt a halo mass 
estimator. We suggest to use the red-band optical lu- 
minosity together with the morphology-dependent mass- 
to-light ratios to estimate the relative mass of galaxies. 
An example of using r-band luminosity as the mass es- 
timator for dark halo plus galaxy systems, can be found 
in Park & Choi (2009), where the method turned out to 
work quite well in the sense that galaxy properties show 
interesting dependence on local and global environments 
at physically meaningful scales. Note that we don't need 
to know the absolute value of halo mass if parameters 
are normalized by their RMS values as we do here. Us- 
ing any halo mass estimator monotonically proportional 
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Fig. 6. — The x-component of peculiar velocity estimated from 
halo number density (top panel) and from halo mass density (bot- 
tom) compared with that from matter density. 

to the actual halo mass will improve the accuracy of the 
resulting potential field. 

3.3. Redshift- space Distortion 

When galaxy distance is obtained from redshift, the 
galaxy distribution becomes biased in such a way that 
clusters and groups are stretched, filaments appear more 
prominent by broadening the interior but compressing 
the exterior, and voids look elongated along the line of 
sight (Kaiser 1987). These redshift space distortion ef- 
fects generate error when the gravitational potential is 
estimated directly from the galaxy distribution in red- 
shift space. We will show here that this error can be 
almost entirely removed by making a linear correction to 
the redshift space distribution of galaxies. 

We calculate the peculiar velocities of dark halos from 
the dark halo number density or mass density in redshift 
space using Equation (11). They are compared with the 
true peculiar velocities in Figure 6. The upper panel 
shows the x-component of the peculiar velocity at each 
galaxy location calculated from halo number density, and 
the bottom one shows that from halo mass density. It 
can be seen that the relation is a little tighter in the 
second case. 

The top panel of Figure 7 shows the correlation be- 
tween the shear magnitudes estimated from halo number 

density and from matter density. The correlation be- 
comes stronger when the halo mass density is used for 
the reconstruction as shown in the middle panel. The 
redshift space distortion effects become tiny when the 
redshift space positions of dark halos are corrected for 
the peculiar velocity, which is calculated from the red- 
shift space distribution of dark halos. This is shown in 
the bottom panel. Ideally one can iterate this correc- 
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Fig. 7. — (top) Magnitude of shear tensor calculated from redshift 
space halo number density versus that from matter distribution, 
(middle) Same, but for redshift space halo mass density, (bottom) 
Same, but for redshift space halo mass density with the linear 
correction for the redshift effects. The numbers at the lower right 
corner are correlation coefficients. 

tion process, but we found that one time correction was 
sufficient. 

3.4. Environmental Dependence of Halo Properties 

In previous sections we studied how accurately one 
can estimate the gravitational potential and its functions 
when only biased and finite samples are available. Now 
we will examine how various physical properties of dark 
matter halos depend on environment including gravita- 
tional potential or its functions. The purpose is to see if 
their dependence on gravitational potential or its func- 
tions is different from that on local density. So it is nec- 
essary to see the correlation between gravitational po- 
tential and local density at galaxy positions. Figure 8 
shows the potential (top), ellipticity (middle), and pro- 
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Fig. 8. — Correlation of local density with gravitational potential 
(top), ellipticity (middle), and prolateness (bottom) of shear ten- 
sor. All are estimated from the dark matter density field smoothed 
over 6?t-'^Mpc. 



lateness of the shear tensor (bottom panel) as a function 
of overdensity. Each variable is normalized by its stan- 
dard deviation. It demonstrates that, even though there 
exists an overall correlation with the local density, the 
dispersion of the potential or shear tensor at fixed local 
density is very large. This reflects the fact that the grav- 
itational potential picks up large scale power (see Eq. 1), 
and thus has information on environment independent of 
local density. 

In Figure 9 we shows four physical parameters of dark 
matter halos as a function of local density, potential, el- 
lipticity, and prolateness of the shear tensor. Points are 
the dark halos with mass greater than 2.9 x lO^^/i~^M0. 
Dashed and solid lines are the mean relations for dark 
halos with Mh > 2.9 x 10" and 1 x IO^^/i^^Mq, respec- 
tively. Only 1/800 of dark halos in the simulation are 
plotted, but the mean relations are obtained from 1 /8 of 
dark halos. As expected, massive halos form in high den- 
sity regions (top-left panel) because of high interaction 
and merger rates (Park, Gott, & Choi 2008; Fakhouri & 
Ma 2009). This plot very much resembles the local den- 
sity dependence of galaxy luminosity obtained from the 
CfA survey (Fig. 12 of Park et al. 1994) and the SDSS 
survey (Fig. 8 of Park et al. 2007). It also tells that 
the maximum mass which a dark halo can acquire is a 
function of the local density. In low density regions mas- 
sive halos simply cannot form by the present epoch. A 
similar observation can be made from the panel showing 
the velocity dispersion of dark halos (the left panel of the 
third row). The mean halo mass is relatively higher in 
low potential regions or in high shear ellipticity and pro- 
lateness regions, but its dependence on these parameters 
is weaker than that on local density. This dependence 
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Fig. 9. — Relations between physical properties of dark halos and environmental parameters for all types of halos with mass above 
Mfi > 2.9 X 10^^h~^ Mq. Considered halo properties are mass (top row), spin parameter (second row), velocity dispersion (third row), 
and short-to-long axis ratio (bottom row) . The left column shows the relation as a function of overdensity in the dark matter density field 
smoothed over 6/i~^Mpc. In the second, third and fourth columns gravitational potential, cUipticity and prolateness of the shear tensor 
are used as environmental parameters, respectively. Dots are halos with mass M/^ > 2.9 X 10^^ Mq, and the solid lines are the mean 
relations. Dashed lines are for halos with Mf, > 1.0 X lO^^h'^ Mq. 



of halo mass on potential and shear is partly due to the 
correlation of potential and shear with local density. 

The second row of Figure 9 shows dependence of the 
spin parameter A on the environmental parameters. The 
spin parameter is defined by (Peebles 1969, Gardner 
2001, and for an alternative definition see. Bullock et 
al. 2001, Shaw et al. 2006, Hetznecker & Burkert 2006) 

A = Jvirl^^virr/'/GMvif, (12) 

where Mvir, Jvir and E'vir sltc the mass, total angular mo- 
mentum and energy of a dark halo, respectively. Con- 
trary to other studies (Maccio et al. 2007) we detect 



dependence of A on all four environmental parameters. 
For halos of Mh > IO^^H^^Mq the spin parameter is 
about 0.034 in very low density regions, but increases 
to about 0.039 in high density regions. Dark halos also 
tend to have higher spin in high shear magnitude regions. 
However, it should be noted that dark halo spin is not 
sensitive to any of these environmental parameters. In 
particular, dependence of the spin on the tidal shear is 
weak, showing only about 4% difference as the environ- 
ment changes from low to high shear ellipticity regions 
(solid line in the third panel of the second row of Figure 
9). 




Fig. 10. — The spin parameter A and axis ratio c/a for isolated halos (top two rows), central halos (middle two rows), and satellite halos 
(bottom two rows) as functions of environmental parameters used in Figure 9. Dots in each panel are halos with mass 3 X 10^^h~^yiQ < 
< 5 X lO^'^/i~^M0 and the long-dashed lines are their mean relations. The short-dashed lines are for halos with 3 X 10^^ h~^M.Q < 
Mh <5 X lO"fe~lM0, and solid lines are for halos with 3 X lO^^/i-^M© < Mh < 5 X lO^^h-^MQ. 
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Fig. 11. — Dependence of c/a on the halo mass [left) and on the spin parameter (right). The filled circles, filled boxes, open circles, and 
the open triangles show the mean distributions of the FoF, central, isolated, and satellite halos. Gray dots show the scatter plots of 30,000 
halos selected in the FoF sample. 
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The panels at the bottom of Figure 9 show that dark 
halos are more spherical in high density, low gravitational 
potential, or high shear regions. Such a trend is stronger 
for more massive halos. On average, high-mass halos are 
preferentially more spherical than low-mass ones, and the 
difference is larger in high density /shear regions. This is 
consistent with the observational finding of Park et al. 
(2007) that early-type galaxies in the SDSS survey are 
rounder as their luminosity increases (see their Figure 
15). However, it contrasts with an opposite finding in 
simulations by many authors (Avila-Reese et al. 2005; 
AUgood ct al. 2006; Maccio et al. 2007, 2008) who re- 
ported that low-mass halos are more spherical than high- 
mass halos at a given epoch in their simulations. The 
main reason for this discrepancy is the difference in the 
halo definition as shown below. 

To study the dependence of halo properties on envi- 
ronmental parameters in a more sensible way we divide 
our halo sample into subsets according to halo mass and 
type. We consider three halo types: isolated, central, and 
satellite halos. Since halos can have a different growth 
history depending on their type, it will be interesting 
to see the halo property-environment relation separately. 
We are also going to fix the halo mass to subtract the 
halo mass dependence of a parameter from its environ- 
mental dependence. Figure 10 shows the spin param- 
eter and axis ratio for isolated (top two rows), central 
(middle two rows), and satellite (bottom two rows) ha- 
los. Dots in the top two rows arc the isolated halos with 
3 X lOi2/i-iM0 < Mh <5x lOi2/i-iM0 and the long- 
dashed lines are the mean relations. Only 1/25 of halos 
are plotted. The short-dashed lines are for the halos 
with 3 X lO"/i-iM0 < M,, < 5 X IO^/j-IMq, and solid 
lines are for the halos with 3 x 10^^/i~^Mq < Mh < 

5 x lO^^/i~^M0. We plot only the spin and c/a param- 
eters here because the velocity dispersion shows almost 
no dependence on environment in these subsamples with 
fixed halo mass ranges, where c/a is the ratio of the mi- 
nor and major radii in a triaxial shape (Chandrasekhar 
1969; de Zeeuw & Franx 1991; Jing & Suto 2002; Smith 

6 Watts 2005; AUgood et al. 2006). In the middle and 
the bottom two rows of Figure 10 the points and lines 
correspond to the same mass ranges as in the top two 
rows but they are for central and satellite halos, respec- 
tively. 

A few observations can be made from Figure 10. First, 

the dependence of A and c/a on the potential and shear is 
still very weak (sec the second, third, and fourth columns 
of Figure 10). When the parameters are studied as a 
function of local density, we find diverse relations that 
depend not only on halo mass but also on halo type. An 
interesting finding is that the spin parameter of satellite 
halos tend to be constant or even decrease as the local 
density increases while that of isolated or central halos 
increases. It seems that the spin of satellites has de- 
creased as they interact more frequently with other halos 
including the central halo in higher density environment. 
An observational evidence for this interpretation is that 
the spin of late-type galaxies decreases as they approach 
their neighbors within the virial radius (Cervantes-Sodi 
et al. 2010). There is a clear trend for the spin parameter 
to decrease as the halo mass increases for all halo types. 
The spin and local density tend to be more positively 
correlated for more massive halos. 



It is also found in Figure 10 that halos are more spher- 
ical in higher density regions for all three types. Further- 
more, the dependence of c/a on local density is stronger 
for more massive halos. Following the same interaction 
picture, we interpret this phenomenon as a result of tidal 
interactions between neighboring halos. Halos in high 
density regions will suffer from the strong tidal force of 
neighboring halos more frequently, and can become more 
spherical on average. Our results are consistent with 
Maccio et al. (2007) who found isolated and central ha- 
los are more spherical in high density regions when the 
density is smoothed by an 8/i^^Mpc tophat even though 
their main conclusion was independence of halo proper- 
ties on large-scale background density. 

The dependence of halo shape on halo mass is not 
monotonic. Figure 10 shows halos are more spherical as 
mass increases from ~ 4 x lO-'^M© to 4 x lO^^M©, then 
more elongated as mass increases to ~ 4 x lO^'^M©. This 
result might seem dissonant with some previous reports 
that more massive halos are less spherical (Avila-Reese 
et al. 2005; Allgood et al. 2006; Maccio et al. 2007). The 
disagreement may originate from difference in halo type 
and mass range. Avila-Reese et al. (2005) and Maccio et 
al. (2007) used the FoF or spherical overdensity (SO) al- 
gorithm to identify halos, respectively. Therefore, a halo 
in these works can be a group of halos connected or em- 
bedded with one another and can have miiltiple density 
peaks. The massive halos located at the intersections 
of dense filaments are more likely to be connected with 
other halos along the filaments by tenuous bridges of par- 
ticles and appear more elongated. On the other hand, the 
halos in our work (the PSB halos) have one and the only 
one density maximum, and arc separated into isolated, 
central, and satellite halos depending on their positions 
and mass. The PSB halos are to be directly identified 
with galaxies individually (Kim et al. 2008). Therefore, 
the shape of the PSB halos docs not necessarily the same 
as the FoF or SO halos as the halo shape depends on halo 
definition (Bett et al. 2007). 

The correspondence of our result with the previous 
ones can be checked by measuring the shape of the FoF 
halos that are identified in our halo finding before the 
step identifying the PSB halos. Figure 11 shows the re- 
lations between c/a and halo mass, and those bctwcien A 
and c/a for FoF (dots and filled circles), isolated (open 
circles), central (boxes), and satellite (open triangles) ha- 
los. The FoF halos are most spherical when their mass 
is near 2 x lO^^/i^^M©. Note that c/a has no maximum 
for isolated halos and is increasing as the halo mass in- 
creases. The spin is a monotonically decreasing function 
of c/a for all types of halo. Figure 11 agrees very well 
with Figure 12 and 13 of Bett et al. (2007). 

The overall message of Figure 10 seems that the in- 
ternal physical properties of dark halos arc mainly con- 
trolled by small-scale physics. This is supported by the 
fact that halo properties depend most sensitively on the 
local density rather than on the gravitational potential. 
It is also consistent with the finding of Park & Choi 
(2009) and Park & Hwang (2009) that galaxy properties 
in the general field and clusters are significantly affected 
by gravitational and hydrodynamic interactions with the 
nearest neighbor galaxy. However, as Figure 8 shows, 
the gravitational potential has information independent 
of that of the local density, and it is worth to explore if 
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the parameters like the gravitational potential or local 
shear can tell something more about galaxy formation. 



4. SUMMARY AND DISCUSSION 

In this paper we demonstrated the gravitational po- 
tential and its functions can be reasonably accurately 
estimated from an observational sample that covers only 
a finite volume of the universe and suffers from tracer 
biasing and redshift space distortion. We found that the 
error in the gravitational potential and shear tensor de- 
creases rapidly as one moves inside the survey bound- 
aries. In the case of shear tensor the error becomes less 
than 20% of its RMS value in the regions separated from 
the survey boundaries by more than about SO/i^^Mpc. 
This requires the sample size to be much larger than 
60/i~^Mpc for an environment study with accurate es- 
timation of the potential and its functions. Our study 
also shows that the effects of halo biasing on the gravi- 
tational potential estimation can be greatly reduced by 
weighting dark halos (or galaxies) by their mass as was 
done by Park et al. (2008) and Park & Choi (2009) in 
their studies of small and large-scale density dependence 
of galaxy properties. Accuracy in the estimation starts 
to fall down rapidly when the halo mass cut is larger 
than lQ^^h~^MQi. This means that the mass density 
and potential fields estimated from the distribution of 
the Luminous Red Galaxies (LRGs) will be quite inac- 
curate because the mean separation of a volume-limited 
sample of the SDSS LRGs is about 20/i^iMpc (Gott et 
al. 2009), which corresponds to the halo mass cut of 
about 2 X W^SH-^Mq. 

The error due to the redshift space distortion effects 
can be also reduced by using dark halo mass density in 
the estimation of potential. But even more reduction 
can be achieved by correcting the observed (i.e. redshift 
space) distribution of dark halos for the peculiar veloc- 
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ity. It was sufficient to use the peculiar velocity linearly 
estimated from the redshift space distribution of dark 
halos. After making the peculiar velocity correction and 
using halo mass weight the error due to the redshift space 
distortion becomes tiny. 

We showed there exists large dispersion in the gravita- 
tional potential and the shear at fixed local density. It 
demonstrates the potential has large-scale information 
independent of local density. We inspected the depen- 
dence of dark matter halo properties on local density, 
gravitational potential, shear ellipticity and prolateness. 
Among these environmental parameters the local density 
shows the strongest correlation with the internal physi- 
cal parameters of dark halos. When halo mass is fixed, 
the spin and shape parameters are nearly independent of 
the potential and shear tensor but depend sensitively on 
the background density in the case of massive halos, in 
particular. 

In the following paper we will analyze the Main Galaxy 
sample of the SDSS DR7 catalog to examine the depen- 
dence of various galaxy properties on these environmen- 
tal parameters. This will extend our understanding on 
the environmental effects on galaxy formation and evo- 
lution. 



The authors thank the anonymous referee for valu- 
able comments. HBP thanks Prof. Hyung-Mok Lee for 
useful comments. JHK appreciates the support of the 
College of Applied Science, Kyung Hee University re- 
search professor fellowship in CAS-RP-KJH2009. CBP 
acknowledges the support of the Korea Science and Engi- 
neering Foundation (KOSEF) through the Astrophysical 
Research Center for the Structure and Evolution of the 
Cosmos (ARCSEC). 



Hwang, H. S. & Park, C. 2009, ApJ, 700, 791 

Hockney, R. W. & Eastwood, J. W. 1981, Computer Simulations 

Using Particles, McGraw-Hill, New York 
Hogg, D. W., et al. 2004, ApJ, 601, L29 
Hubble, E. & Humason, M. L. 1931, ApJ, 74. 43 
Jing, Y. P. k Suto, Y. 2002, ApJ, 574, 538 
Jing, Y. P., Suto, Y., k Mo, H. J. 2007, ApJ, 657, 664 
Kaiser, N. 1987, MNRAS, 227, 1 
Kauffmann, G., et al. 2004, MNRAS, 353, 713 
Kim, J. & Park, C. 2006, ApJ, 639, 600 
Kim, J., Park, C, & Choi, Y.-Y. 2008, ApJ, 683, 123 
Kuehn, F. & Ryden, B. S. 2005, ApJ, 634, 1032 
Lee, J., & Lee, B. 2008, ApJ, 688, 78 

Lee, J. H., Lee, M. G., Park, C, & Choi, Y.-Y. 2010, MNRAS in 
press 

Lee, J. & Pen, U.-L. 2000, ApJ, 532, L5 

Lee, J. & Pen, U.-L. 2002. ApJ, 567, Llll 

Lemson, G. & Kauffmann, G. 1999, MNRAS, 302, 111 

Maccio, A. V., Dutton, A. A., van den Bosch, F. C, & Frank, C. 

2008, MNRAS, 391, 1940 
Maccio, A. V., Dutton, A. A., van den Bosch, F. C, Moore, B., 

Potter, D., & Stadel, J. 2007, MNRAS, 378, 55 
Monaco, P. & Efstathiou, G. 1999, MNRAS, 308, 763 
Navarro, J. F., Abadi, M. G., & Steinmetz, M. 2004, ApJ, 613, L41 
Oemler, A. 1974, ApJ, 194, 1 
Park, C. & Kim, Y.-R. 2009, arXiv:0905.2268 
Park, C. & Choi, Y.-Y. 2009, ApJ, 691, 1828 
Park, C, Gott, J. R., & Choi, Y.-Y. 2008, ApJ, 674, 784 
Park, C, Choi, Y.-Y., Vogeley, M. S., Gott, J. R., & Blanton, M. R. 

2007, ApJ, 658, 898 
Park, C, Vogeley, M. S., Geller, M. J., & Huchra, J. P. 1994, ApJ, 

431, 569 

Park, C, Gott, R. J., & da Costa, L. N. 1992, 392, L51 
Peebles, P. J. E. 1969, ApJ, 155, 393 

Porciani, C, Dekel, a., & Hoffman, Y. 2002, MNRAS, 332, 325 



12 



Park, Kim & Park 



Postman, M., & Geller, M. 1984, ApJ, 281, 95 
Reid, B. &c Spergel, D. N. 2009, ApJ, 698, 143 

Rojas, R. R., Vogeley, M. S., Hoyle, F., & Brinkmann, J. 2004, 
ApJ, 617, 50 

Shaw, L. D., Weller, J., Ostriker, J. P., & Bode, P. 2006, ApJ, 645, 
815 

Sheth, R. K., et al. 2003, ApJ, 594, 225 

Skibba et al. 2009, MNRAS, 399, 966 

Smith, R. E. & Watts, P. I. R. 2005, MNRAS, 360, 203 



Spergel, D., et al. 2007, ApJS, 170, 377 

Vitvitska, M., Klypin, A. A., Kravtsov, A. V., Wechsler, R. H., 

Primack, J. R., & Bullock, J. 2002, ApJ, 581, 799 
White, S. D. M. 1984, ApJ, 286, 38 
de Zeeuw, T. & Franx, M. 1991, ARA&A, 29, 239 
Zheng, Z., Zehavi, I.. Eisenstein. D. J., Weinberg, D. H., & Jing, 
Y. 2008. larXiv:0809. 1868 



